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Summary. — The problem of construction of a continuous (macroscopic) matter 
model for a given point-like (microscopic) matter distribution in general relativity 
is formulated. The existing approaches are briefly reviewed and a physical analogy 
with the similar problem in classical macroscopic electrodynamics is pointed out. 
The procedure by Szekeres in the linearized general relativity on Minkowski back- 
ground to construct a tensor of gravitational quadruple polarization by applying 
Kaufman's method of molecular moments for derivation of the polarization tensor 
in macroscopic electrodynamics and to derive an averaged field operator by utiliz- 
ing an analogy between the linearized Bianchi identities and Maxwell equations, is 
analyzed. It is shown that the procedure has some inconsistencies, in particular, (1) 
it has only provided the terms linear in perturbations for the averaged field oper- 
ator which do not contribute into the dynamics of the averaged field, and (2) the 
analogy between electromagnetism and gravitation does break upon averaging. A 
macroscopic gravity approach in the perturbation theory up to the second order on 
a particular background space-time taken to be a smooth weak gravitational field 
is applied to write down a system of macroscopic field equations: Isaacson's equa- 
tions with a source incorporating the quadruple gravitational polarization tensor, 
Isaacson's energy-momentum tensor of gravitational waves and energy-momentum 
tensor of gravitational molecules and corresponding equations of motion. A suitable 
set of material relations which relate all the tensors is proposed. 
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1. Introduction 

In using the Einstein equations for a matter distribution in the form of a set of 
point-like mass constituents, there is a problem of adequate application, or validity, of 
the Einstein equations when such a matter distribution is substituted by a continuous 
matter distribution while the field operator in the left-hand side of the equations is kept 
unchanged. This problem as it stands in cosmology( 1 ) is called the averaging problem [2, 
3, 4, 5, 6]. Indeed, let us consider the Einstein equations in the mixed form( 2 ) 

(1) «Tr £/3 - l^fTr^ = -< discrctc) 
with 

(2) ^ (discrctc) (,) =5>>Z[* - Z (i )(T (i) )] 

i 

where is a energy-momentum tensor for a point-like mass moving along its world line 
z M = (r(j) ) parameterized by T(j) and i counts for the matter particles in the distribu- 
tion (2). Changing the discrete matter distribution to a continuous (hydrodynamical) one 
in the right-hand side of (1), which is the standard approach in cosmology [2, 3, 4, 5, 6] 
made phenomenologically on the basis of assumption about the uniformity and isotrop- 
icity of distribution of galaxies, or cluster of galaxies, throughout the whole Universe, 
means an implicit averaging denoted here by (•) 

(3) tf^\x) - Tf**°\x) = (£t w ° fi [x *0)(t«))]^ • 

Given a covariant averaging procedure (•) for tensors on space-time, the averaging out 
of (1) with taking into account (3) brings 

(4) (9 a ^) - \^r, v ) = -KTf Y ^°) . 

An important point regarding the averaged equations (4) is that in this form they are 
just algebraic relations between components of the smoothed hydrodynamical energy- 
momentum tensor and the average products of the metric tensor by the Ricci tensor 
(g aer ep) an d cannot therefore be taken as field equations. By splitting the products out 
as (g ae r<:p) = (g ae )(r e p) + Cp where is a correlation tensor, the averaged equations 

(4) become 

(5) (g ae )(r, ) - ^ W><»» - - K r; (hydro) - c% + \^c\. 

( 1 ) For discussion on the other physical settings on general relativity facing the same problem 
see [1]. 

( 2 ) The mixed form is preferable here for the reason that it contains only products of metric 
by curvature. On contrary, the covariant or contravariant forms of the Einstein equations have 
triple products of metric by metric by curvature. 
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Here {g a ') and (r a( g) denote the averaged inverse metric and the Ricci tensors which 
are supposed to describe the gravitational field due to the matter distribution j^( h y dro ) 
A simple important observation [7, 8] now is that the averaged Einstein equations (5) 
are still not "real" field equations - just a definition of the correlation tensor as a 
difference between (4) and (5). The origin of this fundamental fact is that the average 
of the non-linear operator of (1) on the metric tensor g pG is not equal in general( 3 ) to an 
operator of the same form on the average metric (g pa )- 

(6) ^ [g at r^ - ^<rr„„) [ ffp .]) ± (V><^) - ^ W><r>>) [{g P «)\ ■ 

In order to return them the status of the field equations one must define the object 
and find its properties using information outside the Einstein equations. 

To resolve the averaging problem, and to consider it in a broader context as the 
problem of macroscopic description of gravitation, the approach of macroscopic gravity 
has been proposed [5, 7, 8, 11, 12, 13, 14, 15] (see [7, 8, 13] for discussion of the problem 
and references therein, [1] for discussion of the physical status of general relativity as 
either a microscopic or macroscopic theory of gravity). A covariant space-time volume 
averaging procedure for tensor fields [5, 11, 14], has been defined and proved to exist 
on arbitrary Riemannian space-times with well-defined properties of the averages. Upon 
utilizing the averaging scheme, the macroscopic gravity approach has shown that (i) 
averaging out Cartan's structure equations brings about the structure equations for the 
averaged (macroscopic) non-Riemannian geometry and the definition and the properties 
of the correlation tensor Cp, (ii) the averaged Einstein's equations (5) become then 
the macroscopic field equations and they must be supplemented by a set of differential 
equations for the correlation tensor, (Hi) it is always possible to extract the field operator 
of the form (1) for the Riemannian macroscopic metric tensor G pl/ and its Ricci tensor 
M pv with all other non-Riemannian correlation terms going to the right-hand side of 
(5) to give geometric correction to the averaged energy-momentum tensor j i ^( h y dlo \ 
It is been also shown [7, 8, 15] that only in case of neglecting all correlations of the 
gravitational field the averaged equations (5) becomes the macroscopic Einstein equations 
with a continuous matter distribution 

(7) G at M cP - \^G^M^ = -«T; (hydro) , 

which reveals the physical status of using the standard procedure in cosmology of claiming 
(5) to be the Einstein equations (7) after substitution of the matter model (3). The 
physical meaning, dynamical role and magnitude of the gravitational correlations must 
be elucidated in various physical settings. There is some evidence that they cannot be 
negligible for cosmological evolution (see, for example [16] for an estimation of the age 
of Universe in a second order perturbation approach). 

( 3 ) It should noted that the inequality (6) has been observed in all possible averaging settings, 
for example, for a volume space-time averaging in [2, 5], in the framework of a kinetic approach 
in [9] and for a statistical ensemble averaging in [10]. Relations between different averaging 
procedures are discussed in [7, 8]. 
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2. Macroscopic media in general relativity 

Derivation of the macroscopic (averaged) Maxwell field operator in macroscopic elec- 
trodynamics is easily accomplished due its linear field structure and the main problem 
consists in the construction of models of macroscopic electromagnetic media (for example, 
diamagnetics, magnetics, waveguides, etc.) [17, 18, 19, 20], which relates to the structure 
of the averaged current. In general relativity the problem of construction of macro- 
scopic gravitating medium models is hardly elaborated due to the following reasons: (a) 
existing mathematical and physical difficulties in establishing the form of the averaged 
(macroscopic) operator in (5) for the field equations of macroscopic gravity recedes the 
interest in development of macroscopic gravitating media; (b) posing on its own more or 
less realistic problem with a discrete matter creates mathematical and physical problems 
due to nonlinearity and non-trivial geometry of gravitation, to mention, for example, the 
iV-body problem, the problem of statistical description of gravity, etc. and (c) being re- 
lied on physically motivated phenomenological arguments (uniformity isotropy, staticity, 
etc.) most applications of general relativity deal with effective continuous media if even 
a starting physical model is discrete in its nature like in cosmology (see Section 1) or in 
description of extended bodies in general relativity (see [1] for discussion of the physical 
status of general relativity). 

The kinetic approach in physics is known to provide a general scheme for introduction 
of characteristics of continuous media with a known distribution function of a discrete 
configuration. But the advantages of such generality are often greatly weakened in partic- 
ular applications by difficulties of solving the Boltzmann equation to find a distribution 
function of interest. This applies to a great extent to general relativity where despite the 
formulation of the general relativistic Boltzmann equation [21, 22] the kinetic approach 
still remains useful for general definitions and considerations rather than being a working 
tool (see, for example, [9]) to derive a specific model of a macroscopic medium. 

In case of the macroscopic electrodynamics together with the volume space-time av- 
eraging on Minkowski space-time the formalism of statistical distribution functions has 
been utilized (see [19] and references therein) and it is of importance for the mathe- 
matically well-posed derivation of the macroscopic theory and the general structure of 
averaged current starting from microscopic electrodynamics of point-like moving charges. 
Further application of the macroscopic theory requires usually mainly phenomenologi- 
cal considerations to establish material relations between macroscopic average fields and 
induction field necessary to make an overdctermined system of macroscopic equations 
determined. A correct derivation of material relations is known to require [23] averaging 
the microscopic equations with a given microscopic matter model during accomplishing 
an averaging of the microscopic field equations. Though it is the only self-consistent way, 
the elaboration of such kind of approach still remains a challenge even for simple physical 
settings. 

On the other hand, volume (space, time, space-time) averaging procedures maintain 
their importance, direct physical meaning in application to macroscopic settings and 
their extreme clearness and descriptiveness. A volume averaging is also known to be 
unavoidable in all macroscopic settings (including statistical approaches) [20, 23, 24, 25] 
and space-time averages of physical fields are known to have the physical meaning as 
directly measurable quantities [26, 27] (for discussion see [7], [8] and references therein). 
That greatly motivates and supports interest in applying approaches with various av- 
eraging schemes in physics despite corresponding (mostly mathematical, not physical) 
difficulties in the rigorous formulation of averaging procedures. 



GRAVITATING MACROSCOPIC MEDIA IN GENERAL RELATIVITY AND MACROSCOPIC GRAVITY5 



The paper aims to approach the problem of construction of gravitating macroscopic 
media in general relativity by using an appropriate space-time averaging scheme. 

3. Szekeres' gravitational polarization tensor 

The approach of Szekeres [28] to construct a tensor of gravitational quadruple po- 
larization is based on the following assumptions: (a) the linearized theory of gravity 
on Minkowski background; (b) the linearized field equations are taken as the linearized 
Bianchi identity to employ an analogy between gravitation and electromagnetism; (c) 
the covariant method of molecular moments of Kaufman is applied to construct a tensor 
of quadruple gravitational polarization. 

The equations under consideration are the contracted Bianchi identities 

(8) C^upty' kJ^^p 

where C^ upa is the Weyl tensor interpreted as free gravitational field [29, 30], J^ vp is a 
kind of "matter current" for the energy-momentum tensor t^ u 

(9) Jpvp — J[pu]p = — {tp[fj,;u] ~ p[pt ,v\) ■> 

(10) <W P = o. 

equations (8) are analogous the Maxwell equations 

(11) fpS = ^Jp 

with (10) being comparable with the conservation of the electromagnetic current 

(12) V = 0. 

Let us consider a number of particles labeled by i and having masses nrii and which 
are moving in there own effective gravitational field along world lines (rj). A physical 
parameter which characterizes such a distribution is a typical characteristic distance I 
between neighbouring particles. Then the corresponding microscopic energy- momentum 
tensor has the form 

dTi dn 

Assume now that due to gravitation the particles form into groups, a kind of gravitational 
molecules, which will be labeled by index a. From the physical point of view that means 
the presence of another parameter which is a characteristic size (diameter) L of such 
a molecule with L I. It is a longwave macroscopic parameter and its presence in a 
microscopic system will be defining the dynamics of the system on the distances of order 
of L. The microscopic energy-momentum tensor (13) becomes now 

iUa f rW Jy» Jy v 
(14) t^^{x)=C- X Y J Y. J m ^J^fV - y^r a ) - S t (T a )]dT a 
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where y^(r a ) is a world line of the a-th molecule center of mass [28] and s^(r a ) = 
z i( T i) ~ y^{ T a) is a vector connecting the i-th particle with the center of mass of the 
molecule including this particle. Let us apply now the method of molecular moments of 
Kaufman [31] to represent (14) as a series expansion in powers of sf under assumption 
that the effective gravitational field which is created by moving gravitational molecules 
is a weak field and the perturbations of the gravitational field due to relative oscillations 
of gravitating particles in molecules are small compared with the mean effective field. 
After averaging out over a typical size of the gravitational molecule (for an averaging 
procedure see [5, 11]) one gets( 4 ) in accordance with the Szekeres procedure 

(15) ^(molec)^ = T (frcc)^ + D ^P p + QW ^ 

where T^ ee ^ v is the energy-momentum tensor of molecules, which has the form similar 
to (13) with substitution i by a, D^ vp is the tensor of gravitational dipole polarization 
that can be incorporated into the quadrupole term, and Q^ P<T the tensor of gravitational 
quadrupole polarization 



(16) 



Q » vp o = c -i { J2 f q ^5\x - y a )dr a ), 



which has the symmetries of the Riemann tensor. The expression for the covariant 
gravitational quadrupole moment q% upa is defined as 

(17) 
where 

(is) <r = £™^?<, 



2 dr a (ds\ ds P v 



< 19 > ^-fS'-^-V*..- dr." 

(20) K ~ 3 ^ mi dn \ dr a Si dr a Si dn Si S * dr a ) ' 

Upon averaging (8) over the typical size of gravitational molecule the following equations 
were obtained: 

(21) <cv^r - «(^ cro) > 

where 

/99\ i T(micro)\ _ /j-( micro )\ _i_ 1 „ /-/-(micro) \ 

(22) {J^p ')--{t p[ll )M + oVpIuP M> 



( 4 ) No explicit averaging procedure had been used in [28] and averaged relations and equations 
were being written rather on the basis of heuristic considerations than a rigorous analysis. 
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or 

(23) Ppvpa = ~^{~Qpcre[p' ~ ~^ r ]p[p.Q~' rrje' ) ,v] 

and 

(OA) I 7(micro)\ _ T-(free) p ,o 

\") \ ,J pvp I ~ ° jxvp r pvpa ■ 

The expression (24) is analogous to the expression for the averaged electromagnetic 
current (j( mlcro )^} for a bunch of charged particles moving along their world lines in 
the effective electromagnetic field in accordance with the microscopic equation (11) with 
j (micro) n wnen particles are grouped into molecules [31] 

(25) (j( micr °)^ _ j(hee)p, _ c ppv ^ 

where the polarization tensor P^ v is defined as an average of the quadruple polarization 
moment of the molecules (see [31, 28] for details) 

(26) ^ = <£ f dr a p^8\x-y a )). 

a J 

Then equations (21) can be rewritten as the macroscopic equations( 5 ) 

(27) E^ vpa ' = nJp Up 
for the gravitational induction tensor E^ vpa defined as 

(28) E^upa — {Gpupo) ~t~ ^Ppupo- 

The macroscopic equations (27) are analogous to the Maxwell macroscopic equations 
obtainable by means of averaging the microscopic equations (11) with j( mlcro )M w ith 
taking into account (25) 

4-7T 

(29) H" v v = — j( frcc )^ 

c 

for the electromagnetic induction tensor defined as 

(30) = (fn+^Pp,- 

Unfortunately, at this point the analogy between the electromagnetism and gravitation 
which holds on the level of (8), (10) and (11), (12) breaks. Indeed, the formal similarity 
of (27), (28) and (29), (30) does not possess the structural analogy between averaged 
electromagnetism and gravitation: (A) the gravitational induction tensor E^ vpcr does 
not have any more the symmetries of the Weyl tensor compared with keeping the 
symmetries of f^; (B) it is constructed from the second derivatives of the polarization 



( ) Gravitational macroscopic equations similar to (21) are known to have been proposed first 
in [32]. 
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tensor Qp Vpa compared with the linear algebraic structure of the electromagnetic induc- 
tion tensor H plJ in terms of the polarization tensor P^ v - it is thus impossible to proceed 
with the formulation of phcnomenological material relations between Ep Vpa and (C pvpa ) 
as possible in electromagnetism (relations between H pu and (/ M „), or amongst the fields 
E, D, B, H and J). 

Even more important issue is that analysis of the macroscopic field equation (21) 

(31) {Cfjtupu) P = K(Jp U p = Kj^vp ^ P^vpa ,<7 > 

O(e) O(l) 0(e 2 ) 

where e is a parameter measuring the value of deviation from the flat space, requires 
one to put into agreement the orders of magnitude of all quantities and reveals that the 
linearized Weyl tensor should be zero under averaging (C pi/pa ) = 0. The Weyl tensor 
must be estimated in the perturbation theory up to the second order as it was done for 
the matter current jjw P CI °^ in the left-hand side of (8). So considering some physical 
features of the polarization tensor neither the macroscopic equations (21), nor any other 
field equations had in fact been employed [28]. 
On the basis of the expression 

(32) (i£r r °)> = T^r ] + \q wv « p ° 

the following material relations have been suggested 

(33) Qmjo = (Gij)N = egdojo 

where N is the average number of molecules per unit volume, dj is the quadrupole 
moment of a molecule 

(34) dj = J p(x)SxiSxjd 3 x, 

p = p(x) is the matter density in molecules, Sxi is a vector between neighbouring particles 
of a molecule. The quantity e g has been called the gravitational dielectric constant and 
in Newtonian approximation found to be 



1 mA c 

35 e g - —N 



where A is the average linear dimension of a typical molecule, m is the average mass of 
the molecules, ujq is a typical frequency of harmonically oscillating particles in molecules. 

4. Macroscopic gravity equations 

The gravitational field created by the particles is defined by Einstein's equation 

(36) <Tr e/J -^%, = -/<^ 
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The Einstein equations (36) for the distribution (14) are of the form: 

(37) <Tr e/3 - ^£<rr> = -< molec) . 

Averaging the left-hand side of the Einstein equations (37) following the Isaacson's 
high-frequency approximation approach [33, 34], using the averaging procedure [5, 11] 
(one can also use the Isaacson's averaging procedure [33, 34], see also [12]) and with 
taking into account (16) brings the averaged Einstein equations in the form: 

(38) <) - \ 9 $RW = -«(T£-> + T(«w) + l -Q^- P °l 

where T^ W ^ is Isaacson's energy-momentum tensor of gravitational waves [33, 34]. All 
members in equation (38) can be shown to be of the same order of magnitude 0(1/ L 2 ). 
The macroscopic equations (38) give the equations of motion for molecules( 6 ) 

(39) T (free)»v v = 

conservation of the energy-momentum of gravitational waves 

(40) T (GW) "% = 0, 
and an identity for the gravitational polarization 

(41) Q^pa'"^ = 0. 

The system of equations (38)-(41) is underdetermincd - there are 20 unknown compo- 
nents of the tensor of gravitational polarization. It is possible to formulate two natural 
material relations. The first relation is between the traceless part of the quadrupole 
polarization tensor 



(42) Qp,pva — Qppva .Qp,vPpa-> 



1 

4'' 

where P pa = Q^p^a, and traceless energy-momentum tensor of gravitational waves 



(43) \Qpp^ = \T^\ 

where A = \(x). Relation (43) can be shown to be always valid in the geometrical optics 
limit. 

The second material relation connects the remaining part of the polarization tensor 
Qppua, its trace P prJ , with a projection of the curvature tensor on the world line of an 
observer (electric part of the curvature tensor) 



(44) Ppa e -R-\ipua 



( 6 ) Under assumption that the background metric in the left-hand side of (38) represents a weak 
gravitational field on the flat background one can use the covariant derivatives with respect the 
metric in all relations instead of partial derivatives with respect to the flat metric. 
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where u M is the observer 4-velocity (4- velocity of the molecule centre of mass) and e = 
e(x). The relation (44) can be shown to lead to the correct expression for the 3-tensor 
of the average quadrupole gravitational moment (34) so that 

(45) = (P Q0 = 0, P 0l = 0, Pij = (Gij)N) . 
Then the material relation (33) can be recovered in the form 

(46) Qiojo = (Gij)N — e g R i0 jo 

that gives e — e g with the gravitational dielectric constant e g defined by (35). 

Thus the system of equations (38)-(40), (43), (44) is fully determined and can be used 
to find the gravitational and polarization fields for the macroscopic gravitating systems. 

* * * 

Roustam Zalalctdinov would like to thank Remo Ruffini for hospitality in ICRA where 
the work has been done in part. 
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